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ABSTRACT 

We present radio, optical, ultraviolet, and X-ray observations of the isolated, thermally-emitting neu- 
tron star RX J0720.4— 3125 using the Parkes radio telescope, the Very Large Array, the Hubble Space 
Telescope, and the Chandra X-ray Observatory. From these data we show that the optical/UV spectrum 
of RX J0720.4— 3125 is not well fit by a Rayleigh- Jeans tail as previously thought, but is instead best fit 
by either a single non-thermal power-law or a combination of a Rayleigh- Jeans tail and a non-thermal 
power-law. Taken together with the X-ray spectrum, we find the best model for RX J0720.4— 3125 
to be two blackbodies plus a power-law, with the cool blackbody implying a radius of 11-13 km at 
an assumed distance of 300 pc. This is similar to many middle aged (10 5-6 yr) radio pulsars such as 
PSR B0656+14, evidence supporting the hypothesis that RX J0720.4— 3125 is likely to be an off-beam 
radio pulsar. The radio data limit the flux at 1.4 GHz to be < 0.24 mJy, or a luminosity limit of 
4ird 2 F < 3 x 10 25 

^300 er S s s 1 > an( i we see n0 s i§ n °f extended nebulosity, consistent with expectations 
for a pulsar like RX J0720.4-3125. 

Subject headings: pulsars: individual (RX J0720.4— 3125) — stars: neutron — X-rays: stars 
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1. INTRODUCTION 

Thermally-emitting neutron stars have been the targets 
of many observations recently, as these sources can poten- 
tially reveal the equation of state (EOS) of neutron stars, 
and thereby explore nuclear physics in realms inaccessible 
from laboratories (Lattimer & Prakash 2000). To obtain 
the EOS from the spectrum seems simple: determine the 
effective angular size from spectral fits, multiply by the 
distance (obtained from other means), and one has the 
apparent radius. This radius can be converted into the 
physical radius through use of mass. The radius is the cru- 
cial quantity in differentiating between EOS, as most EOS 
predict a distinctive but small range of radii for a large 
range of masses. However, in order to use a neutron star 
to determine the EOS, one needs to (1) be certain that the 
radiation is from the surface, and (2) have a thorough un- 
derstanding of this emission. Radio pulsars and accreting 
binaries arc unsuitable since non-thermal magnctosphcric 
emission or emission from the accreting material far ex- 
ceeds the surface emission. 

Therefore the identification of the nearest neutron stars 
by ROSAT was a major advance in the field (see reviews 
by Motch 2001 and Treves et al. 2000). Most of these 
sources do not have significant non-thermal emission, so 
they are prime targets for studies leading to the EOS. The 
closest of these sources, RX J1856.5— 3754, has been the 
subject of much inquiry lately for just this purpose (e.g., 
Walter & Lattimer 2002; Drake et al. 2002; Braje & Ro- 
mani 2002). 

RX J0720.4-3125 was discovered by Haberl et al. (1997) 
as a soft (kT ~ 80 eV), bright X-ray source in the ROSAT 
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All-Sky Survey — the second brightest neutron star that 
is not a radio pulsar. Given its very low hydrogen col- 
umn density (N H ~ 1 x 10 20 cm" 2 ), nearly sinusoidal 
8.39-s pulsations, relatively constant X-ray flux, and very 
faint (B = 26.6 mag), blue optical counterpart (Kulka- 
rni & van Kerkwijk 1998; Motch & Haberl 1998), it was 
classified as a nearby, isolated, thermally-emitting neu- 
tron star. It is perhaps the second closest source (next 
to RX J1856.5— 3754) that does not show significant non- 
thermal emission. While originally thought to be an ac- 
creting source or possible an old magnetar, recent tim- 
ing measurements limit the original magnetic field to be 
smaller than 10 14 G, eliminating the magnetar hypothe- 
sis. However, the discoveries of radio pulsars with peri- 
ods longer than 4 s (Camilo et al. 2000; Young, Manch- 
ester, & Johnston 1999) have led to the suggestion that 
RX J0720.4— 3125 is instead an off-beam radio pulsar, 
likely with a magnetic field at the high end of the radio- 
pulsar range (B ~ 10 13 G; Kaplan et al. 2002b; Zane et al. 
2002). 

The spectrum emerging from a thermally-emitting neu- 
tron star depends significantly on the composition of the 
surface (e.g., Romani 1987). In the past, three models have 
generally been considered. The first, a blackbody, is sim- 
ple but not physically motivated. Next, light element (H 
or He, possibly due to accretion from the ISM) atmosphere 
have few features, all in the (extreme) ultraviolet, and 
peak, for a given temperature, at a substantially higher 
energy than a black-body. Finally, heavy element (Fe, Ni) 
atmospheres peak at a similar location to blackbodies but 
have many spectral features at a variety of wavelengths. 

None of these models can reproduce the X-ray and op- 
tical data for RX J0720.4-3125 and RX J1856. 5-3754 
(Pons et al. 2002; Walter & Lattimer 2002; Drake et al. 
2002; Paerels et al. 2001). The X-ray data are well fit by 
blackbodies, but these blackbodies underpredict the opti- 
cal flux. The H/Hc models also match the general shape, 
but they overpredict the optical flux and the implied radii 
are larger than is possible for a neutron star. The Fe mod- 
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els have too many lines and edges to match the X-ray 
spectra. Consequently, the current generation of single- 
component models have been unsuccessful in fitting both 
the X-ray and optical fluxes while having radii consistent 
with those of neutron stars (Pons et al. 2002; Kaplan, van 
Kerkwijk, & Anderson 2002c). We see, though, how both 
X-ray and optical data are required to fully constrain the 
models. 

Motivated thus, we present observations of RX J0720.4— 3125 
from radio to X-ray wavelengths aimed at determining its 
spectral energy distribution and from that its underlying 
properties (composition, magnetic field, and geometry). 
In Section 2 we present new optical/UV data from the 
Hubble Space Telescope, and undertake detailed modeling 
of the optical/UV spectrum. In Section 3 we present spec- 
troscopic data from the Chandra X-ray Observatory, and 
in Section 4 we present searches for radio sources (both 
persistent and pulsating) with the Very Large Array and 
the Parkes radio telescope. Finally, in Section 5 we dis- 
cuss the spectrum of RX J0720.4— 3125, and present our 
conclusions in Section 6. In the following, all radii re- 
fer to the radiation radius as observed at infinity, 
(i?oo = i? P hys/v/ 1 ~~ 2GM/i? phys c 2 , where i? phys is the 
physical radius). 

2. OPTICAL AND UV DATA 

2.1. Observations and Analysis 

We observed RX J0720.4-3125 with the Space Tele- 
scope Imaging Spectrometer (STIS) aboard the Hubble 
Space Telescope (HST) four times, covering wavelengths 
from 125 nm to 900 nm; the observations are summarized 
in Table 1. 

2.1.1. Optical Data 

The optical data (50CCD mode) consist of eight un- 
filtered CCD observations taken in a four-point dither 
pattern. We assembled the images using the drizzle al- 
gorithm (Fruchtcr & Hook 2002) giving a plate-scale of 
25.3 mas pixel -1 . We show the stacked image in Figure 1. 

We performed standard aperture photometry on the 
stacked STIS image using IRAF's daophot package. The 
sky level was estimated using an annulus from (X'75-l'.'OO. 
The source flux was measured within an aperture of radius 
0'.'5. While there were no aperture corrections strictly ap- 
propriate for a source with the color of RX J0720.4— 3125 
(B — V i=a —0.3 mag), we used the bluest of the color- 
dependent aperture corrections available (T. Brown 2002, 
personal communication) to correct the flux to an infinite 
aperture. We estimate that the aperture correction intro- 
duces an uncertainty of < 0.02 mag, as blue sources like 
RX J0720.4— 3125 have less scattered light than redder 
sources and therefore the aperture corrections are better 
determined (see Kaplan et al. 2002a) . Given the very wide 
50CCD bandpass (FWHM m 441 nm) a single zero-point 
flux is not appropriate for all source spectra. Therefore, as 
a first order estimate, we calculated the zero-point flux at 
the mean wavelength of the filter (given in Table 1), assum- 
ing an input spectrum with Fx oc A~ 4 and Ay ~ 0.1 mag, 
appropriate for RX J0720. 4-3125 (Kulkarni & van Kerk- 
wijk 1998). 

For astrometric purposes we used 10-s and 60-min R- 
band images taken with the Low-Resolution Imaging Spec- 




FlG. 1— STIS optical image of RX J0720.4-3125. North is 
up, East to the left. The scale bar in the upper right indicates 5". 
RX J0720.4— 3125 is indicated by the arrow, and sources A, B, and 
C from Haberl et al. (1997) are also marked. 



trograph (Okc et al. 1995) on the 10-m Keck I telescope 
(the observations are described in Kulkarni & van Kerk- 
wijk 1998). We determined the centroids of 231 stars from 
the Guide Star Catalog version 2.2 (GSC-2.2) on the 10-s 
image, rejecting 50 objects that were overexposed or ap- 
peared to be incorrectly identified. The pixel coordinates 
were corrected for instrumental distortion using a cubic 
radial distortion function 5 , and we then fit for the zero- 
point position, rotation, and plate-scale of the image. The 
rms was 0'.'12 in each coordinate. We then fit the 70 non- 
saturated stars from the 10-s image to the 60-min com- 
posite image, again correcting for instrumental distortion. 
This fit had rms of 0'.'024 in each coordinate. Finally, we 
performed a fit using 25 stars from the 60-min image to 
determine the zero-point, scale, and rotation of the STIS 
image (which had been corrected for distortion during the 
drizzle process), giving an rms of 0'.'025 in each coordinate. 
Overall, the STIS image is tied to the frame of the GSC-2.2 
with uncertainty of O'.'Ol in each coordinate, or tied to the 
International Coordinate Reference Frame (ICRF) with 
uncertainty of about 072 in each coordinate. The final po- 
sition for RX J0720.4-3125 is (J2000, epoch MJD 52106) 
a = 07 h 20 m 24!961, 5 = -31°25'50'/21. This supersedes 
the position of Kulkarni & van Kerkwijk (1998), as the 
significance of the new detection is far higher and the po- 
sition is directly tied to the ICRF. The HST position is 
consistent with the X-ray position which has uncertainties 
of w 076 in each coordinate (Kaplan et al. 2002b). 

2.1.2. UV Data 

For the UV MAMA data, we corrected the arrival times 
of the data to the Solar System barycenter using the 
stsdas task odelaytime. We then filtered the photon lists 
for the FUV data for periods of high background. These 
occurred at the beginning of each orbit as HST was go- 
ing into the Earth's shadow. Using a background annulus 
from 3'.'0-3'.'5, we estimated the median background level 
from a light-curve binned to 50-s and only used the data 

5 See http://alamoana.keck.hawaii.edu/inst/lris/coordinatcs.html. 



The Spectrum of RX J0720.4-3125 
Table 1 

Summary of STIS Observations 
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2001-Jul-16 


5342 


5148 


1.56 


0.102 




26.68 ±0.10 


NT IV MAMA /V^Otv 


9001 Alio- fl^ 


OOU\J 


2286 


9 ^8 


97(S 


n 09Q 


oq so _i_ n i 4 


FUV MAMA/25SrF 2 


2002-Jan-28 


4800 


1447 


2.80 


0.340 


0.087 


21.97 ±0.11 


FUV MAMA/25MAMA 


2002-Feb-13 


3850 


1360 


2.97 


0.353 


0.078 


21.6 ± 0.2 










Calculated for 


a„ = 1.40 


and Aq = 0. 


10 mag. f 


CCD/50CCD 






5367 


1.41 


0.102 




26.68 ±0.10 


NUV MAMA/25Qtz 






2310 


2.56 


0.270 


0.029 


23.83 ±0.14 


FUV MAMA/25SrF 2 






1450 


2.79 


0.339 


0.088 


21.96 ±0.11 


FUV MAMA/25MAMA 






1365 


2.95 


0.352 


0.079 


21.7 ±0.2 



a Effcctive wavelength and normalized extinction; see Appendix A of van Kerkwijk & Kulkarni (2001b). 

b Aperture correction from a radius of Ci.'5 to infinity. See § 2 for details. 

c Change in the zero-point magnitude due to degradation of the MAMA detectors. 

d Magnitude in the STMAG system, with msT = —21.1 — 2.51og 10 i 7 A, corrected to infinite aperture. 

° Appropriate for the RJ fit in Table 4. 

'Appropriate for all fits except the RJ fit in Table 4. 



Note. — All filter curves were taken from the synphot database. Ao = 4500 A. 



where the background level was within ±3 a of the median. 
We found that the background was noticeably high only 
for the FUV MAMA/25SrF 2 data, where we have retained 
4800 s of the original 5500 s. 

We then performed aperture photometry on the raw 
photon data. The source flux was computed using a radius 
of 075, where the signal-to-noise was relatively high and 
the aperture corrections were well-defined. The raw source 
and background counts are given in Table 2. The STIS 
MAMAs have wide point-spread- functions (psfs), with 
substantial flux beyond 0'.'5, and therefore aperture correc- 
tions are particularly important. We took the monochro- 
matic aperture corrections for the MAMAs (T. Brown 
2002, personal communication) convolved with the ex- 
pected source spectrum and the filter throughputs to com- 
pute aperture corrections for each filter, which we give in 
Table 1. Another issue with the STIS MAMAs is that the 
sensitivity changes with time (Stys et al. 2002) at the level 
of a few percent per year. Using a wavelength-dependent 
fit to the sensitivity changes (D. Stys 2002, personal com- 
munication; R. Diaz-Miller 2002, personal communica- 
tion; these corrections are now incorporated into the STIS 
pipeline), we computed weighted zero-point corrections for 
the NUV and FUV MAMA data that are listed in Table 1. 
These corrections were small, less than 0.1 mag. As with 
the optical data, the zero-point fluxes were calculated at 
the mean wavelengths of the filters. 

The STIS MAMA data arc time-tagged with 125 ^s res- 
olution, but we could not detect the 8.39-s periodicity 
present in the X-rays. We find 90% confidence upper lim- 
its on the Fourier power (normalized to have unity mean 



and rms) present at the X-ray period (given by Kaplan 
et al. 2002b) of 5.3 and 5.9, for the NUV and FUV/SrF 2 
data sets, respectively. These translate into limits on the 
rms pulsed fraction of 0.16 and 0.19, respectively. 

2.2. Spectral Fits 
2.2 A. Power-law Fits 

We already know from Kulkarni & van Kerkwijk (1998) 
that the optical spectrum of RX J0720.4— 3125 is roughly 
approximated by the Rayleigh- Jeans portion of a black- 
body curve: F\ cx A~ 4 . We can now investigate this 
more quantitatively and over a wider range of wavelengths. 
Here we use the optical and UV data presented in this 
paper as well as the B and R photometry from Kulka- 
rni & van Kerkwijk (1998); for the data from Kulkarni & 
van Kerkwijk (1998), we used the zero-point fluxes from 
Bessell, Castelli, & Plez (1998). 

Our spectral fitting followed van Kerkwijk & Kulkarni 
(2001b). We fit the data with a sum of extincted power- 
law (PL) of the form 

-(2+a„) 

1Q -A A' x /2.5^ (1) 



F X =F 



where av is the spectral index 6 , Fq is the observed flux 
at the reference wavelength Ao, Aq is the extinction at 
Ao, and A' x is the normalized extinction at wavelength A 
(A' x = A\/Ao). We use the reddening curve of Cardelli, 
Clayton, & Mathis (1989), with corrections to the optical 
and UV portions from O'Donnell (1994). We chose A = 

6 From this definition wc also have the more standard F v oc u a " . 
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Table 2 

Source and Background Counts for STIS UV Data 



Detector/Filter 


Counts 


Source 


Background 


NUV MAMA/25Qtz 


3534 ± 59 


2788 ± 15 


FUV MAMA/25SrF 2 


611 ±25 


128 ±3 


FUV MAMA /25MAMA 


14502 ± 120 


13788 ±58 



Note. — No aperture corrections have been applied. 
Background counts are normalized to the same area as 
the source counts. 



4500 A (Ao = 1.29 Ay, or 4 - 4), as this is the mean 
wavelength of the data. 

For the fit, we calculated likelihood values as we varied 
Fq, Ao, and a v . However, as the spectral shape changes 
(i.e. variations in Ao and a v ) the aperture corrections, 
reference wavelengths, extinctions, and zero-point correc- 
tions (from Table 1) also change (these changes are most 
significant for the STIS/50CCD data). Therefore, dur- 
ing the iterations of the fitting, we recomputed all of the 
spectrum-dependent quantities for each combination of A 
and a v . 

The results of the fitting for the two basic models - 
Rayleigh- Jeans 7 (RJ) and unconstrained power-law (PL) 
are given in Table 3. These fits give values of the ex- 
tinction Ao ~ 0.5 mag. This is much higher than ex- 
pected from other observations. From the X-ray spec- 
trum of RX J0720.4-3125, we know that the hydro- 
gen column density is Njj ~ 1-3 x 10 20 cm~ 2 , which 
implies Ay « 0.07 mag (Predehl & Schmitt 1995) or 
Ao ~ 0.09 mag. Kulkarni & van Kcrkwijk (1998) have 
placed an upper limit on the reddening of Eb-v < 0.04, as 
RX J0720.4— 3125 is in the foreground of the open cluster 
Collinder 140, which implies Ao < 0.10 (for the standard 
ratio of Ry = 3.2). We therefore expect small values of 
the extinction: Ao < 0.15 mag, allowing for uncertainties 
in Nh, Rv, and the relation between Nh and Ay. From 
this we can reject the models in Table 3. 

Therefore we constrained Ao from the information 
above. To formally include this in our fit, we performed a 
maximum likelihood fit with the following prior distribu- 
tion for Ao: 

f Ao (A ) = A/" 6 (A)|0.09 mag, 0.06 mag), (2) 

where Mn(x\fj,, a) is a generalized Gaussian distribution of 
degree n (n is even) with the form: 

Af n (x\n,a) = — " — r " ' 



exp 



(- 

Vn 



V2cr 



(3) 



For n = 6, this distribution essentially requires that Ao 
be between mag and 0.18 mag. We could have used a 

7 For blackbodies of the temperatures considered here, the Rayleigh- 
Jeans approximation holds at the shortest UV wavelength used to 
better than 4% — considerably smaller than the measurement error. 



uniform prior with Ao = 0-0.15 mag, but the sharp edges 
of this distribution can make for discontinuities in the re- 
sulting posterior distributions, so we opted for Equation 2, 
which is a smoothed version of the uniform distribution. 

We consider four power-law models for the fit. The first, 
given in the PL column in Table 4, is a single power-law 
fit to the optical/UV data alone. The second, given in 
the RJ column in Table 4, is a single power-law fit to 
the optical/UV data, but where the power-law index is 
that of a Rayleigh- Jeans tail (a v = 2). The third is a 
fit where there are two power-laws, given in the PL±RJ 
column in Table 4: one has an unconstrained index, and 
the other has a v — 2. Finally, the fourth fit, given in the 
PL±X-ray column in Table 4, has a Rayleigh-Jeans power- 
law present, but its normalization is set by the X-ray fit 
(§ 3; the uncertainties in the X-ray flux extrapolated to 
optical/UV wavelengths is w 10%). 

In Tables 3 and 4 we give % 2 values for each fit. These 
values were computed from the data without taking into 
account the prior distributions used in Table 4. The best- 
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Fig. 2.— Optical/UV spectrum of RX J0720.4-3125, with data 
from Kulkarni & van Kerkwijk (1998) and this paper. The PL fit 
is shown, with: the best-fit model (dashed line), the best-fit model 
corrected for extinction (dotted line) and the extrapolation of the 
Chandra blackbody fit (dash-dotted line; "CXO BB"; see § 3). The 
different bands arc labeled. 
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Table 3 

Fits to RX J0720. 4-3125 Optical/UV Data with A Unconstrained 

Parameter Type of Fit 

PL" PJ 



A (mag) 0.47(15) 0.59(6) 



1.9(2) 2" 



F (xl0~ 19 ergs s" 1 cm" 2 A -1 ) 2.5(6) 3.0(3) 
Bp (mag) b 26.0(2) 25.82(10) 

X 2 o n 

DOF 3 4 

X 2 /DOF 0.3 0.3 

a The spectral index such that F x oc A~ (2+Qv) . 
b B-band Vega- magnitude corresponding to Fq. 



Note. — Numbers in parentheses are 68% confidence lim- 
its in the last digit(s). Values without confidence limits 
were held fixed for the fit. Values with subscript are for 
A = 4500 A. 



Table 4 

Fits to RX J0720.4— 3125 Optical/UV Data with A Constrained 

Parameter Type of Fit 

PL RJ PL+RJ PL+X-ray 



Ap (mag) 0.10(6) 0.19(2) 0.10(4) 0.09(4) 





„ a 


2 A- 1 ) 


1.40(4) 


2 


0.3(7) 


1.14(12) 


F (xl(T 


19 — 1 - 

ergs s cm 


1.50(12) 


1.12(7) 


0.8(2) 


1.24(9) 




Bo (mag) b 




26.58(8) 


26.88(7) 


27.4(2) 


26.79(8) 












2 


2 


F ,2 (xlO 


-19 —1 

ergs s cm 


- 2 A" 1 ) 






0.6(2) 


0.27 




Bo.a (mag) b 








27.5(3) 


28.5 




x 2c 




1.9 


26.5 


0.8 


1.4 




DOF 




3 


4 


2 


3 




X 2 /DOF 




0.6 


6.6 


0.4 


0.5 



a The spectral index such that F\ cx A _(2+a " ) . 
b B-band Vega-magnitude corresponding to Fo- 

c The x 2 values are raw values that do not take into account the prior distributions 
(e.g., Eqn. 2). They are there only as a guide, showing which models do and do not 
fit the data independent of the prior distributions. 



Note. — Numbers in parentheses are 68% confidence limits in the last digit(s). 
Values without confidence limits were held fixed for the fit. Values with subscript 
arc for A = 4500 A. A was constrained by use of Equation 2. 
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MA) 

Fig. 3.— Optical/UV spectrum of RX J0720.4-3125, with data 
from Kulkarni & van Kerkwijk (1998) and this paper. The PL+RJ 
fit (unconstrained Raylcigh-Jeans tail plus a second PL) is shown: 
the best-fit model (dashed line), the best-fit model corrected for 
extinction (dotted line), the two components of the model (solid 
lines), and the extrapolation of the Chandra blackbody fit (dash- 
dotted line; "CXO BB"; see § 3). The different bands are labeled. 



fit values of the parameters were computed not by using 
the x 2 values themselves buy through the marginalized 
likelihood functions that incorporated the prior. The \ 2 
values are there only as a reference, to show that even 
with the prior distribution in place the fits are still good. 
One should not use the typical A\ 2 technique (Press et al. 
1992, p. 697) for determining parameter uncertainties - 
while the PL models from Tables 3 and 4 have x 2 's that 
differ by 0.9, the PL model from Table 3 is excluded very 
significantly by the fit in Table 4, as shown by the small un- 
certainties in the parameters of Table 4. Instead the con- 
fidence limits given in Tables 3 and 4 are single-parameter 
68% limits determined from the marginalized likelihood 
functions. 

We find a good fit for the single unconstrained power- 
law, shown in Figure 2. This model has x 2 = 1.9 
for 3 degrees-of- freedom (DOF). The Rayleigh- Jeans fit, 
though, is significantly worse, with x 2 = 26.5 for 4 de- 
grees of freedom under the same prior assumptions. We 
can therefore reject a v — 2 with > 98% confidence. For 
the PL+RJ fit, the results are given in the PL+RJ col- 
umn of Table 4. The fit, shown in Figure 3, is good, with 
X 2 = 0.8 for 2 DOF. For the final fit, PL+X-ray, we also 
find an acceptable value of \ 2 — 1-4 for 3 DOF. We can 
conclude that all of the fits except the RJ fit are accept- 
able. 

As mentioned above, the values of (A), A' x , the aperture 
correction and the zero-point correction change depending 
on the spectral model. However, except for the RJ fit, all 
of the models have sufficiently similar flux distributions, 
despite the different contributions from different compo- 
nents, that the values are essentially the same for these 
models. We give the values of these parameters for the 
best-fit values of a v and A in Table 1; these will apply 
to the PL, PL+RJ, and PL+X-ray models. The only fil- 
ter whose calibration changes significantly is the extremely 
broad-band STIS/50CCD. 



As seen in Figure 2, the PL fit has a shallower slope 
than the extrapolation of the Chandra blackbody spec- 
trum. While an extrapolation of the optical/UV PL does 
not intersect the X-ray spectrum, it does come to within 
a factor of 1.2 (at 142 A). By the lower energy end of 
the Chandra data (0.1 kcV w 125 A), the optical/UV PL 
must have turned over as it is not seen in the X-rays. The 
power-law component of the PL+X-ray fit behaves simi- 
larly. 

The best-fit Rayleigh- Jeans component of the PL+RJ 
model is a factor of 2.4(4) above the extrapolation of the 
X-ray blackbody (Fig. 3 and Tab. 4). The two components 
contribute equally at A = 4930 A: at shorter wavelengths 
the Rayleigh- Jeans component dominates, while at longer 
wavelengths the non-thermal component dominates. The 
non-thermal PL component is above the X-ray extrapola- 
tion in the optical regime, intersects it at 2680 A, and the 
continues below it. Therefore, the non-thermal PL would 
not be seen in soft X-rays or in the radio (Fig. 4) , although 
it may approach the X-ray spectrum at energies > 2 keV. 

2.2.2. Disk Fits 

We also considered fits to the optical/UV data that in- 
clude a disk of accreting material, such as that proposed 
to account for the X-ray luminosity and periodicity of 
RX J0720.4-3125 (e.g., Wang 1997; Koncnkov & Popov 
1997; Alpar 2001). For the disk spectrum we used the 
model of Perna, Hcrnquist, & Narayan (2000). We found 
fits using a disk model to be unsatisfactory. As there are 
too many free parameters to do a formal fit (the inner 
and outer disk radii, the disk inclination, as well as two 
undetermined efficiency factors), we varied subsets of the 

Energy [keV) 




■•• (Hz) 



Fig. 4.— Broad-band spectrum of RX J0720.4-3125, from ra- 
dio to X-rays. The absorption-corrected optical/UV data (§ 2) 
are plotted as filled circles, the absorption-corrected X-ray data 
(§ 3) as points, the VLA upper limits (§ 4.1) as filled triangles, 
and the Parkes upper limits for pulse widths of 0.6-3% (§ 4.2) as 
the open triangles. The models arc: X-ray blackbody (thick solid 
line), Rayleigh- Jeans component of the PL+RJ fit to the optical/UV 
data (dashed line), non-thermal component of the PL+RJ fit (dash- 
dotted line), and the overall PL+RJ fit to the optical/UV data (thin 
solid line). The ±1 <r uncertainties on the non-thermal PL are shown 
by the shaded region. 
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parameters by hand. We first considered disks that extend 
in to the corotation radius and out to infinity. For most 
conceivable disk inclinations (i < 85°) the disk alone is a 
factor of 10-20 above the optical data, and this is without 
any contribution to the optical emission from the neutron 
star surface. Toward the short-wavelength end of the disk, 
where the emission decreases below the level of our data 
and the contribution is primarily from the inner edge of 
the disk, the slope is entirely inconsistent with the opti- 
cal/UV data: it goes approximately as F\ oc A 4 , while the 
excess flux in the optical (compared to a Rayleigh-Jeans 
tail) is like F\ oc A~ 2,3 . We therefore consider this disk 
model to be very unlikely for RX J0720.4-3125. 

There are disk models that can reproduce a spectrum 
roughly similar to that observed. This occurs when the 
inner radius is far inside the corotation radius and ap- 
proaches the neutron star surface, while the outer radius 
move inward to ~ 10 8 cm. But while the spectral shapes 
are not inconsistent, the flux predicted by such disks is a 
factor of ~ 100 above the optical/UV data. Also, there 
is no natural reason for the disk to be truncated at such 
small radii (the optical data do not allow for any stellar 
companion). Therefore this disk model is also very un- 
likely for RX J0720.4-3125. 

2.2.3. Variability 

It is possible that the spectrum of RX J0720. 4-3125 
is a Rayleigh-Jeans tail in the optical, and that the devia- 
tions we see are temporal in nature: i.e. RX J0720.4— 3125 
could vary. However, we consider this unlikely First, 
the X-ray flux has been extremely constant over almost 
a decade of observation (Haberl ct al. 1997; Paerels et al. 
2001). Second, similar sources such as RX J1856.5— 3754 
and PSR B0656+14 have exhibited constant optical fluxes, 
again over several years of observations (van Kerkwijk & 
Kulkarni 2001b; Koptsevich et al. 2001). 

Regardless, we can perform a simple test for variabil- 
ity. We have an ongoing series of HST observations that, 
while designed to measure the parallax and proper motion 
of RX J0720.4— 3125, also provide a sensitive flux monitor. 
Only the first two epochs of data have been observed so 
far (at MJDs 52459 and 52532). The data arc from 4950-s 
observations with the Advanced Camera for Surveys on 
HST using the High Resolution Camera (ACS/HRC) in 
the F475W filter. The photometric calibration of ACS 
is not complete, so we cannot directly compare the mea- 
sured flux of RX J0720.4— 3125 to the models presented 
here (although the new data are roughly consistent), but 
we can look for variations in the flux of RX J0720.4— 3125 
and in that of the other sources in the field. We drizzled 
the data using a preliminary model of the ACS/HRC ge- 
ometric distortion, and then performed aperture photom- 
etry on RX J0720.4— 3125 and 11 other sources, ranging 
from much brighter than RX J0720.4— 3125 to about as 
bright as RX J0720.4-3125. The field sources changed 
by 0.02 ± 0.04 mag from the first epoch to the second, 
while RX J0720.4-3125 changed by -0.10 ± 0.09 mag. 
This number is preliminary, but shows that the flux of 
RX J0720.4-3125 changed by at most 10% over two 
months. We will eventually have better-calibrated data 
spanning two years, which will allow us to make a much 
more rigorous test for variability. 



3. X-RAY DATA 

RX J0720.4— 3125 was observed with the Chandra X-ray 
Observatory, using the High Resolution Camera spectro- 
scopic detector (HRC-S) with the Low Energy Transmis- 
sion Grating (LETG); the observations are summarized in 
Table 5. Here we describe the spectroscopic analysis of 
these data — timing analysis is described in Kaplan et al. 
(2002b) and Zane et al. (2002). 

The spectral data were reduced from standard event lists 
using IDL using custom processing scripts as described in 
Marshall & Schulz (2002) and Marshall et al. (2002). Raw 
events were extracted from first and higher orders and 
calibrated with an updated model of the LETGS effec- 
tive area 8 (EA), which was developed from observations 
of PKS 2155-304.. 

Integrating the observed fluxes over the 0.25-3.0 keV 
band gives an observed flux of (9±2) x 10 -12 ergs cm -2 s _1 
and an absorbed luminosity of (9± 2) x 10 31 c?| o er S s s_1 > 
where d = 300d 300 pc is the distance to RX J0720.4-3125 
(Kaplan et al. 2002c). The data were rebinned adaptively 
to provide a signal-to- noise ratio (S/N) of 5 in each bin 
over the 0.10-2.0 keV range. The spectrum, shown in 
Figure 5, was first estimated using the first-order EA only. 
The contributions to the observed counts from high orders 
are estimated by folding a model for first-order through 
the high-order EA (important only below 0.2 keV). 

Following previous analyses (Haberl et al. 1997; Paerels 
et al. 2001), we modeled the continuum by an absorbed 
blackbody; the fitted parameters are given in Table 6. We 
exclude the data below 0.15 keV from the fit where un- 
certainties in the high-order grating efficiencies can be im- 
portant. This fit gives a temperature of 81.4(13) eV and 
a bolomctric luminosity of 2.1 x 10 32 d| 00 ergs s _1 , consis- 
tent with the ROSAT and XMM analyses (Haberl et al. 
1997; Paerels et al. 2001), given the uncertainties in mod- 

This effective area is available 

at http:/ /cxc. harvard.edu/cal/Links/ 

Letg/User/Hrc_QE /EA /correct.ea/letgs.NOGAP_EA_001031 .mod. 




0.1 1.0 

Energy (keV) 



Fig. 5.— LETGS spectrum of RX J0720.4-3125. The bin sizes 
have been varied to provide good signals in each energy bin; the 
uncertainties are about 20% everywhere. Solid line: Model consist- 
ing of one blackbody component. High orders do not contribute 
significantly for E > 0.20 keV, while the data at low energies 
(E < 0.15 keV) are best modeled as the result of the sum of high 
orders. 
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Table 5 

Summary of Chandra HRC-S/LETG Observations 



Date 


Exposure 




Counts 


(UT) 


(ksec) 


Order 


Orders ±1 


2000-Feb-01 


5.4 


929 


671 


2000-Feb-02 


26.3 


4584 


3027 


2000-Feb-04 


6.1 


1119 


687 



Table 6 

One and Two Blackbody Fit to the LETG Spectrum 
Parameter Value 



N H (xlO 20 cm' 2 ) 1.32(14) 1.46(14) 

T hot (xlO b K) 9.45(15) 9.47(15) 

Rhot (km) 6.1(3)d 3 oo 6.1(6)d 3 oo 

T cold (xl0 5 K) ••• 3.72(10) 

RcoM (km) 15(17) 

X 2 163.9 163.1 

DOF 148 146 

X 2 /DOF 1.11 1.12 



Note. — Numbers in parentheses are 68% 
confidence limits in the last digit (s). 




40 

Wavelength (A) 

Fig. 6.— Count spectrum of RX J0720.4-3125 obtained with 
the LETGS. A binning of 0.125 A was used to obtain sufficient 
signals per bin to search for narrow features. Heavy dashed line: 
Expected count spectrum from the single blackbody model shown 
in Fig. 5. Light dotted lines: ±1 a uncertainties about the model. 
The residuals are consistent with statistical fluctuations about the 
model. The sharp edges in the model near the 50-70 A range are 
the result of detector gaps. 



tent with the XMM data (Paerels et al. 2001), although 
now the upper limit for the energy of any narrow features 
is somewhat lower, around 0.2 keV. 



1000 



1 00 ^ 
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cling the effective MI' CMS below 0.2 keV and the higher order 
responses. The fit has a reduced x 2 = 1.11 (for 148 DOF), 
acceptable at the 90% level. Other models, such as those 
with non-thermal power-laws or a second blackbody, did 
not improve the fit. 

The best-fit model is plotted against the fitted data in 
Figure 5. The count spectrum (Fig. 6) was binned at 
0.125 A resolution in order to search for narrow spectral 
features against the continuum model. No significant emis- 
sion or absorption features were found: in Figure 7 we give 
3 a upper limits to the equivalent width of narrow-line fea- 
tures in our data. The Chandra data are generally consis- 



Energy (keV) 

Fig. 7. — Estimates of the 3cr limits that can be placed on any 
emission lines whose FWHMs are comparable to the instrument 
resolution, as a function of energy. The computation uses the model 
fitted to the data (Fig. 5) and the effective areas. We assume that 
candidate features arc only 2 bins (0.25 A) wide. The curves are 
rather smooth, except for locations of chip gaps, so one may derive 
limits on broad features using these curves until the scale of the 
feature becomes comparable to that of the instrument calibration 
uncertainties. Limits on absorption features arc identical when there 
are many counts but are systematically larger at the high and low 
ends of the spectrum, where there are fewer than 25 counts per bin. 
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Fig. 8.— VLA maps of the field around RX J0720.4-3125. Left: 
1.4 GHz map from the 1999 observations. Right: 8.4 GHz map. 
RX J0720.4— 3125 is indicated with the circle at the center (the 
radius of the circle is 10"). 



4. RADIO OBSERVATIONS 

4.1. Synthesis Imaging 

We observed RX J0720.4-3125 with the Very Large Ar- 
ray (VLA) once at 8.4 GHz and twice at 1.4 GHz (summa- 
rized in Table 7) in the standard synthesis-imaging mode. 
All data sets were independently calibrated using AIPS, 
but the two 1.4 GHz observations were combined for imag- 
ing. 

For the 1.4 GHz data, we performed imaging and self- 
calibration in difmap. We iteratively cleaned and self- 
calibrated (phase only) until the gain solution converged. 
Uniform weighting was used, yielding a synthesized beam 
with FWHM w 34". An overall gain adjustment was 
added for one IF, effectively correcting for a non-zero spec- 
tral index across the two IF's. No additional amplitude 
self-calibration was necessary. 

After cleaning, we found rms map noise to be 0.08 mJy, 
a factor of <~ 8 higher than the theoretical thermal noise 
but consistent with confusion (Condon et al. 1998). The 
final image (see Fig. 8) shows a 5.0 mJy point source next 
to the position of RX J0720. 4-3125, but we believe that 
this source is unrelated. For reference, the radio source is 
at J2000 a = 07 h 20 m 28?28(2), S = -31°26'09'/9(3), 46" 
away from the nominal position of the source. No point- 
like or diffuse emission from RX J0720.4— 3125 was found, 
which then gives 3 a upper limits to the flux of a point- 
source of 0.24 mJy and to that of an extended source of 
< 0.43 mJy arcmin~ 2 . 

The imaging of the 8.4 GHz data proceeded simi- 
larly. Cleaning and phase self-calibration were done in 
difmap. Again, no source was found at the position of 
RX J0720.4-3125 (see Fig. 8), giving a 3 a flux limit of 
0.12 mJy. 

4.2. Pulsation Searches 

We observed RX J0720-3125 with the 64-m Parkes ra- 
dio telescope on 11 January 2001 in an attempt to de- 
tect pulsed radio emission. The target position was ob- 
served for 10,800 s at center frequencies of 1374 MHz 
and 644 MHz. At the higher frequency we used the 
center beam of the Parkes multibeam receiver to feed a 
3 MHz x 96-channcl multibeam filterbank (see Lyne et al. 
2000). At 644 MHz the front-end was the Parkes 50-cm 
receiver, and the back-end was a 0.125 MHz x 256-channcl 
filterbank. Both observations employed a 1-ms sample pe- 
riod and one-bit digitization. 

The interstellar dispersion toward RX J0720.4-3125 is 



unknown, but we can estimate it with the latest model of 
Galactic electron density (Cordes & Lazio 2002) which pre- 
dicts a dispersion measure DM = 4 pc cm~ 3 at d = 300 pc 
or DM = 30 pc cm~ 3 at d = 500 pc. We therefore 
take DM = 100 pc cm~ 3 as a conservative upper limit 
to the DM (the search was highly insensitive to DM any- 
way, given the long period of RX J0720.4-3125). Both 
datasets were de-dispersed with dispersion measures up to 
100 pc cm~ 3 and searched for periodicities near the known 
X-ray period (Kaplan et al. 2002b; Zane et al. 2002) using 
standard folding and FFT-based techniques. No pulsar- 
like signals were detected. 

Without a detection, we must estimate the limiting flux 
of a signal. While the search is insensitive to DM, the 
shape of the hypothetical pulse profile strongly affects the 
sensitivity. Assuming a pulse duty cycle of w = 1%, our 
8(J detection limits are 0.2 mJy at 644 MHz and 0.02 mJy 
at 1374 MHz, and they scale approximately as 



Smin, 644 MHz = 2.0 * mJy 

V 1-w 



Smin, 1374MHz = 0.18*/- mJy. (4) 

y i-w 

For signals with high w (i.e. few harmonics) the 1.4 GHz 
sensitivity is comparable to that of the VLA observations 
(§ 4.1), where the large bandwidth of the multibeam sys- 
tem compensates for its smaller area and shorter integra- 
tion. However, for very narrow signals (where many har- 
monics are summed) the periodicity search is a factor of 
<~ 7 deeper than the VLA observations. 

5. DISCUSSION 

5.1. The Spectrum 

We have shown in § 2 that the optical/UV spectrum of 
RX J0720.4— 3125 docs not follow a pure Rayleigh- Jeans 
tail. It is possible that the emission (minus the contribu- 
tion of the Rayleigh- Jeans tail of the X-ray spectrum) is 
entirely non-thermal in origin (like the Crab pulsar, where 
the non-thermal emission entirely overwhelms any ther- 
mal component). However, the slope of the spectrum, 
a v = 1.1 is much steeper than that seen for other pul- 
sars, whole spectral indices range from a v = 0.11 for 
the Crab (Sollcrman et al. 2000) to a„ <~ —1 for other 
sources (Zharikov et al. 2002, and references therein). In 
addition, the low spin-down power of RX J0720.4— 3125 
(E < 2.4 x 10 31 ergs s _1 ; Kaplan et al. 2002b) compared 
to sources like the Crab or even middle-aged radio pul- 
sars means that there is no reason to expect this much 
non-thermal emission from RX J0720.4— 3125. 

Of the six other isolated neutron stars with good op- 
tical/UV data, there are three (all < 10 6 yr old and 
within 500 pc) that show evidence for thermal optical 
emission: the pulsars PSR B0656+14 (Pavlov, Welty, 
& Cordova 1997) and Geminga (Martin, Halpern, & 
Schiminovich 1998) 9 , and the nearby isolated neutron star 
RX J1856. 5-3754 (van Kerkwijk & Kulkarni 2001b). [For 
the 10 4 -yr Vela pulsar, the optical emission (Mignani & 

9 Martin et al. (1998) find the PL+RJ model acceptable for Geminga, 
although a single non-thermal PL is also allowed. In addition, they 
require an absorption feature in the spectrum. 
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Table 7 

Summary of VLA Observations 



Date 


Frequency 


Exposure 


Config. 


Beam Size 


RMS 


(UT) 


(GHz) 


(sec) 




(asec) 


(A*Jy) 


1998-Feb-07 


8.4 


6720 


D^A a 


15" x 5"8 


40 


1999-Feb-18 


1.4 


4410 


DnC 


36" x 32" 


80 


1999-Apr-19 b 


1.4 


7380 


D 







Note. — Observations all had 2 x 50-MHz bandwidths. 

a Data were taken while switching from D configuration to A configuration. 

b Processed with the 1999-Fcb-18 observation. 



Caraveo 2001) is dominated by the non-thermal com- 
ponent, but the X-ray spectrum contains thermal and 
non-thermal contributions (Pavlov et al. 2001).] Of 
these sources, RX J1856. 5-3754 and PSR B0656+14 
have thermal X-ray spectra too, like RX J0720.4— 3125. 
We therefore believe that the optical/UV emission from 
RX J0720.4— 3125 is mostly thermal in nature (i.e. a v = 
2), with the deviation from a Rayleigh- Jeans tail aris- 
ing either due to a multiplicative opacity or to an added 
component. We address each of these models separately. 

5.1.1. One Component Model with Absorption 

It is possible that the emission from RX J0720.4— 3125 
is entirely thermal, but that the underlying Rayleigh- 
Jeans tail is modified by a frequency-dependent absorp- 
tion to give the observed spectrum. An opacity k v oc jA 9 
would give the correct result. However, this model is 
artificial, and is entirely contrary to what is seen with 
RX J1856. 5-3754, PSR B0656+14, and the X-ray spec- 
trum of RX J0720.4-3125 itself (§ 2.2.1). While we cannot 
reject this model based on our data alone, comparison with 
other sources makes it unlikely. 

5.1.2. Two Component Model 

A simpler and more physically motivated model for the 
optical/UV spectrum of RX J0720.4-3125 is that the 
emission is composed of significant thermal emission plus a 
non-thermal PL (the PL+RJ model) similar to the spectra 
of PSR B0656+14 and Geminga. We have plotted bright- 
ness temperatures derived from model optical/UV spec- 
tra for these sources and RX J1856.5— 3754 in Figure 9. 
The sources exhibit some variety in their spectra, ranging 
from purely thermal (RX J1856.5— 3754) to largely non- 
thermal (e.g., PSR B0656+14). While RX J0720.4-3125 
shows less non-thermal emission than PSR B0656+14 or 
Geminga, it is otherwise unremarkable. 

The (extinction corrected) 5-band luminosity of RX J0720. 
is Lb = 6.5 x 10 27 o?| 00 ergs s _1 . If we separate the 
non-thermal component, we find Ls, n on-th = 3.2 x 
10 27 d3 00 ergs s _1 . This compares well with the optical 
luminosities of similarly-aged pulsars given in Zharikov 
et al. (2002), which have 27 < log Lb < 28. As such, 
RX J0720.4— 3125 fits in with the pulsar population de- 
spite its low rotational power (E). We note, though, that 



all of the sources in Zharikov et al. (2002) with r > 10 5 yr 
have very similar values of Lb, despite values of E that 
vary by about 3 orders of magnitude, and that the values 
of Lb show no apparent correlation with E. So it appears 
that is not simply related to E, contrary to what is 
seen in the X-ray regime (Becker & Triimper 1997, see 
below). It may be, as suggested by Zharikov et al. (2002), 
that the efficiency of producing optical emission increases 
with time for r > 10 4 yr so that Lb remains relatively 
constant despite variations in E (i.e. Lb = f]{t)E, where 
di](t)/dt > and dE/dt < 0), or it may just be that the Lb 
is not directly coupled to E. A final possibility is that the 
true distribution of Lb is closely tied to E but that we only 
see the brightest sources in this distribution due to obser- 
vational bias: a source with log Lb = 26.5 at 500 pc would 
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Fig. 9. — Brightness temperature Tg vs. wavelength for for 
-RX25)720.4- 3125 (solid line, this paper), RX J1856.5— 3754 (dashed 
line, van Kcrkwijk & Kulkarni 2001b), PSR B0656+14 (dotted 
line, Pavlov et al. 1997), and Geminga (dot-dashed line, Martin 
et al. 1998). The brightness temperature is defined as Tg(A) = 
F v (X)X 2 Q/2k, and the values were computed from model opti- 
cal/UV spectra assuming a radius of 10 km and distances of 300 pc, 
140 pc, 330 pc (Kaplan et al. 2002c; Walter & Lattimcr 2002), and 
160 pc (Caraveo et al. 1996), respectively. Deviations from a con- 
stant Tb indicate non-thermal emission. 
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have B w 31 mag, below modern detection limits. Simi- 
larly, the very well studied RX J1856.5-3754 at w 140 pc, 
the limit to non-thermal emission is XogLs < 26.5. There- 
fore, any other source with log Lb < 26.5 could not be 
detected, and the true distribution of Lb may be poorly 
represented. 

The slope of the non-thermal PL is a v w 0.3. This 
is consistent (within errors) with the slope of the opti- 
cal PL for the Crab, and is not that far from the slope 
of the spectrum of PSR B0656+14 (a v = -0.45 ± 0.26; 
Koptsevich et al. 2001). The non-thermal spectra of neu- 
tron stars have been observed to follow roughly the same 
slope for over ~ 5 orders of magnitude in energy (e.g., 
Koptsevich ct al. 2001), suggesting a single underlying 
mechanism for the non-thermal optical and X-ray power- 
laws. Here, the non-thermal PL would contribute <~ 1% 
of the flux of the blackbody in the LETG band and would 
therefore be difficult to observe: F(lkeV) = (5l 25 ) x 
1Q-13 cr g S s -i cm -2 At higher energies (> 2 keV) the 
non-thermal PL could contribute substantial X-ray flux 
(although the PL is highly uncertain at these energies; 
see Figure 4), flux that was not detected in the Chandra 
or XMM data. The XMM EPIC-pn data roughly give 
a flux of 2 x 10~ 14 ergs s _1 cm" 2 at 1.5 keV, consistent 
with a blackbody and barely consistent with the extrap- 
olated optical/UV power-law. So, for RX J0720.4-3125 
the non-thermal spectrum is not likely to bridge the X-ray 
and optical regimes. 

Becker & Triimper (1997) have found a rough rela- 
tion between the spin-down power and the non-thermal 
X-ray emission of pulsars: Lx,non-th ~ 10~ 3 E (with 
Lx,non-th in the ROSAT band of 0.1-2.4 keV). Applying 
this to RX J0720.4— 3125 we would expect Lx,non-th 
10 28 ergs s -1 , or Fxnon-th ^ 10~ 15 ergs s™ 1 cm~ 2 . This 
is a factor of > 10 4 less than the observed thermal X- 
ray flux in the same band, and a factor of 10 2 less than 
the extrapolated non-thermal emission. However, there is 
considerable uncertainty in both the relation of Becker & 
Triimper (1997) (see, e.g., Possenti et al. 2002) and in the 
extrapolation of the non-thermal PL to the X-ray regime, 
so the difference between Fx, n on-th predicted from the 
spin-down and that predicted from the optical/UV spec- 
trum may not be significant. 

The thermal component in the optical/UV band is, 
like that of RX J1856.5— 3754, above a simple extrapo- 
lation of the X-ray spectrum (for RX J1856.5— 3754 it ex- 
ceeds the X-rays by a factor of <~ 16; van Kerkwijk & 
Kulkarni 2001b). This could just be a matter of tem- 
perature, though. For the sources PSR B0656+14 and 
PSR B1055— 52, the X-ray spectrum is best-fit by a com- 
bination of a power-law and two blackbody components 
(Greiveldinger et al. 1996; Pavlov, Welty, & Cordova 1997; 
Koptsevich et al. 2001; Marshall & Schulz 2002; Pavlov, 
Zavlin, & Sanwal 2002), where the smaller hot portion 
(presumably the polar cap) contributes the majority of 
flux in the traditional X-ray band but the larger cool por- 
tion (along with the non-thermal emission) gives rise to the 
optical/UV flux (Pavlov et al. 2002). In both cases, the 
blackbody components are all hot enough (> 8 x 10 5 K) to 
appear in the X-ray band, and are therefore well modeled. 

For RX J1856.5— 3754, which is closer in temperature 
to RX J0720.4-3125, Braje & Romani (2002) again ap- 



pealed to a multitcmpcrature surface, but here the second 
blackbody is not observed but only inferred (also see Wal- 
ter & Lattimer 2002). It must be too cool for the X-ray 
band (T co id < 5 x 10 5 K), while still hot enough to ap- 
pear as a power-law in the optical/UV band. While not 
very well constrained, this model gave reasonable results 
for RX J1856.5— 3754, including a weak constraint on the 
radius, and is therefore valuable. We now apply this model 
to RX J0720.4-3125. 

We see from Table 6 that a second (unconstrained) 
blackbody component does not appreciably change the X- 
ray fit. However, we can use the goodness-of-fit (given 
by x 2 ) to constrain the combinations of T co id and i? C oid 
that are allowed. The flux of a Rayleigh- Jeans tail goes 
as F v cx R 2 T. We know that the X-ray blackbody has 
-Rhot = 6.1g?3oo km and Thot = 9.45 x 10 5 K, and since the 
thermal fit to the optical/UV data is a factor of 2.4 above 
the X-ray extrapolation, we find 

T cold = 2.4 x 10 5 2 Kj (5) 

\ "300 / 

where the cold radius i? co id = 15-R C oid,i5 k m has been 
taken to be the radius of the neutron star. Such a cool 
blackbody would not have been seen in the X-ray data 
(as shown by the large uncertainties on i? co id in Tab. 6), 
but we can constrain T co id to be w (3.5-5.0) x 10 5 K, or 
-Rcoid ~ (H-13)d3oo km (at roughly 90% confidence). This 
is similar to the temperature of the cold component found 
for RX J1856.5-3754 (Braje & Romani 2002) — not sur- 
prising since in both cases the cold component was forced 
to give a Rayleigh- Jeans tail in the optical while not giving 
significant contribution in the soft X-rays — and the size 
agrees well with estimates for the radius of a neutron star 
(Lattimer & Prakash 2000). 

Cropper et al. (2001) model the phase-dependent hard- 
ness ratio of RX J0720.4-3125 in XMM data, and find 
that it is consistent with a polar-cap model for a large 
range of cap sizes with angular radii of 10-50°. This 
agrees quite well with our findings (i? co id = 15c?3oo km 
corresponds to an angular radius of w 25°). 

Taken together, these observations show that a RJ+PL 
model for the optical/UV spectrum of RX J0720.4-3125 
makes it entirely consistent with being an off-beam radio 
pulsar, one that likely has a cooler blackbody component 
in the extreme UV. 

5.1.3. Constraints on the Magnetic Field 

We find no significant absorption features in the pulse- 
averaged or pulse-phased spectra over the 0.15-0.80 keV 
band (§ 3). Thus, following Marshall & Schulz (2002) and 
Paerels et al. (2001), we can rule out electron and pro- 
ton cyclotron resonance lines in this range. Paerels et al. 
(2001) already rule out the range 0.03 < B 12 < 0.2 and 
50 < < 200, where B 12 = B/{10 12 G). By extending 
the spectrum down to 0.15 keV, we extend the lower limits 
of the excluded ranges of magnetic fields to B\i — 0.015 
and B12 = 25, though the lower limits could increase if the 
bulk of the emission comes from the equatorial zone, where 
the magnetic field is ~ 50% of the polar value, or if the ab- 
sorbing plasma is far off the neutron surface. We can also 
use the lack of features in the spectrum to exclude hydro- 
gen atmospheres for a range of magnetic field strengths. 
Again following Paerels et al. (2001), we can exclude the 
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range Byi > 15. The excluded ranges of B agree with the 
finding that RX J0720.4— 3125 is not a magnetar (Kaplan 
et al. 2002b), but otherwise the B is consistent with either 
of the models discussed in Kaplan et al. (2002b, see also 
§ 5.2), namely that have B 12 ~ 1-10. 

5.2. Radio Luminosity 

At a distance of 300^300 pc, we limit the 1.4 GHz radio 
luminosity of RX J0720.4-3125 to < 3 x 10 25 d§ 00 ergs s" 1 
(L ra d = iircPF), or following the radio pulsar convention 
L' rad < 0.02 mJy kpc 2 (L' lad = Fd 2 ). This is significantly 
below what is expected of radio pulsars with parameters 
(P and P) similar to those of RX J0720.4-3125. For in- 
stance, the two high-i? radio pulsars discovered by Camilo 
et al. (2000) have i rad - 5 x 10 26 ergs s _1 , while the 8-s 
10 12 -G pulsar PSR J2144-3944 has L lad ~ 10 30 ergs s" 1 
(Young et al. 1999). Similarly, the radio luminosity model 
of Arzoumanian, Chernoff, & Cordes (2002) predicts lumi- 
nosities of io 27 ~ 28 ergs s _1 , depending on the value of P. 
So we can see that RX J0720.4— 3125, if it has any radio 
emission, must be beamed away from the Earth. 

If RX J0720.4-3125 is like PSR J2144-3944, then we 
might expect a similarly narrow radio beam of w ~ 0.6%, 
a beam width that agrees well with the extrapolation 
of Rankin (1993). RX J0720.4-3125 may however be 
more similar to the 10 13 G radio pulsar PSR J1814-1744 
(Camilo et al. 2000), which has a significantly wider beam 
(w i=a 3%). In either case, we can expect that the radio 
beam subtends a small solid angle, making the lack of radio 
emission quite credible. For such beams, our upper lim- 
its to the radio luminosity decrease, as the limit from the 
Parkes data for a source with pulse width of w — 0.6-3% 
at 1.4 GHz is L' rad w 0.002 mJy kpc 2 . Even the imag- 
ing (VLA) limit is quite faint, about a factor of 3 fainter 
than the limit for Gcminga (Seiradakis 1992), and a factor 
of - 30 below that of PSR J0205+6449 in 3C 58 (Camilo 
et al. 2002), but the implied limit for a narrow pulse width 
is far below that of all radio pulsars younger than 10 6 yr 
(Motch 2001). 

With such small beams, the radio-quiet population of 
sources like RX J0720.4— 3125 could potentially be very 
large, up to a hundred times the radio-loud population 
(assuming a sharp cutoff in the radio beam). While likely 
invisible to radio observations, such sources are of course 
bright X-ray emitters, and would be visible to ~ 5 kpc 
in a 30 ksec XMM observation. The total numbers of 
such sources (either radio-loud or -quiet) are small, mak- 
ing statistics uncertain, but there could be cooling radio- 
quiet neutron stars in as much as 1% of XMM observa- 
tions. However, these sources would be all but impossi- 
ble to confirm, as there would be few X-ray photons for 
spectral fitting or pulsation searches, and the optical/UV 
counterpart would be extremely faint. 

6. CONCLUSIONS 

We have shown, through a joint analysis of radio, op- 
tical, and X-ray data, that the spectrum of the isolated 
neutron star RX J0720.4— 3125 cannot be fit by a single 
blackbody model. While statistically we cannot rule out a 
model with a single power-law in the optical/UV domain 
and an X-ray blackbody, from a more general perspective 



we believe that the best-fit model is one with three com- 
ponents: a hot (~9x 10 5 K) blackbody on the polar cap, 
a cool (<~ 4 x 10 5 K) blackbody over the whole surface, and 
a weak non-thermal power-law in the optical/UV. This is 
very similar to the spectra of middle aged radio pulsars 
such as PSR B0656+14 and PSR B1055-52, an observa- 
tion that supports the identification of RX J0720.4— 3125 
with an off-beam radio pulsar. 

RX J0720.4— 3125 appears extremely similar to the very 
nearby RX J1856.5— 3754, perhaps with orientation being 
the only difference between them (Braje & Romani 2002). 
We believe it likely that non-thermal emission and/or pul- 
sations will be detected eventually for RX J1856.5— 3754, 
as suggested by Braje & Romani (2002). Both sources 
seem to have spectra primarily composed of featureless 
blackbodies. If we can develop a full understanding of 
such spectra to properly relate the blackbody radii (such 
as that given here) to the true radii, these sources will be 
ideal targets for the determination of the cquation-of-state. 

In the near future, we will obtain HST astrometry 
allowing us to determine the distance and velocity of 
RX J0720.4— 3125, as well as Ha imaging to search for 
bow-shock nebulae (e.g., van Kerkwijk & Kulkarni 2001a) 
that will place significant constraints on alternate mod- 
els such as accretion. In the off-beam pulsar model, these 
data may limit the luminosity of any particle wind from 
RX J0720.4-3125, and will reduce uncertainty in the EOS 
by determining the conversion from solid angle to radius. 
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